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Abstract. Inhomogeneous chemical evolution models of galaxies which try to reproduce the scatter seen in element- 
to-iron ratios of metal-poor halo stars are heavily dependent on theoretical nucleosynthesis yields of core-collapse 
supernovae (SNe II). Hence inhomogeneous chemical evolution models present themselves as a test for stellar 
nucleosynthesis calculations. Applying such a model to our Galaxy reveals a number of shortcomings of existing 
nucleosynthesis yields. One problem is the predicted scatter in [O/Fe] and [Mg/Fe] which is too large compared 
to the one observed in metal-poor halo stars. This can be either due to the oxygen or magnesium yields or due 
to the iron yields (or both). However, oxygen and magnesium are a-elements that are produced mainly during 
hydrostatic burning and thus are not affected by the theoretical uncertainties afflicting the collapse and explosion 
of a massive star. Stellar iron yields, on the other hand, depend heavily on the choice of the mass-cut between 
ejecta and proto-neutron star and are therefore very uncertain. We present iron yield distributions as function of 
progenitor mass that are consistent with the abundance distribution of metal-poor halo stars and are in agreement 
with observed 56 Ni yields of core-collapse supernovae with known progenitor masses. The iron yields of lower-mass 
SNe II (in the range 10 — 20 Mq) are well constrained by these observations. Present observations, however, do 
not allow us to determine a unique solution for higher-mass SNe. Nevertheless, the main dependence of the stellar 
iron yields as function of progenitor mass can be derived and may be used as a constraint for future core-collapse 
supernova/hypernova models. A prediction of hypernova models is the existence of ultra a-element enhanced 
stars at metallicities [Fe/H] < —2.5, which can be tested by future observations. The results are of importance 
for the earliest stages of galaxy formation when the ISM is dominated by local chemical inhomogeneities and the 
instantaneous mixing approximation is not valid. 
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1. Introduction 

The key to the formation and evolution of the Galaxy 
lies buried in the kinematic properties and the chemical 
composition of its stars. Especially old, metal-poor halo 
stars and globular clusters are ideal tracers of the forma- 
tion process. Although many of the properties of the halo 
component and its substructures have been unveiled, it is 
still not possible to decide whether the Galaxy formed by a 
fast monolithic collapse (Eggen, Lynden-Bell & Sandage, 



1962), by the slower merging and accretion of subgalac- 



tic fragments (Searle & Zinn |1978) or within the context 



of a hybrid picture, comb ining aspects of both scenarios. 
Recently, Chiba & Beers ( 2000 ) made an extensive inves- 
tigation to address this question, concluding that a hybrid 
scenario, where the inner part of the halo formed by a fast, 
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dissipative collapse and the outer halo is made up of the 
remnants of accreted subgalactic fragments, best explains 
the observational data. It also seems to be consistent with 
the theory of galaxy formation based on cold dark matter 
scenarios (see e.g. S teinmetz & M tiller 1995 ; Gnedin 
Klypin et al. |1999|; M oore et al. |l999| ; Pearceet al. 
Bekki & Chiba |2000| Navarro & Steinmetz J2000D . 



1996 



1999 



However, the kinematic structure of the halo alone is 
not sufficient to draw a conclusive picture of the formation 
of the Galaxy. Old, unevolved metal-poor halo stars allow 
us to probe the chemical composition and (in)homogeneity 
of the early interstellar medium (ISM) and its evolution 
with time, since element abundances in the stellar atmo- 
spheres of those stars directly reflect the chemical com- 
position of the material out of which they formed. It is 
almost impossible to determine the age of single stars (ex- 
cept in a few cases where radioactive thorium or uranium 
was detected, see e.g. Cayrel et al. 2001). Therefore, the 
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metallicity Z or iron abundance [Fe/H] of a star is taken 
as an age estimate, knowing that an age-metallicity rela- 
tion can only be used in a statistical sense for the bulk of 
stars (see e.g. Argast et al. 200C, hereafter Paper I). 

Common chemical evolution models mostly assume 
that the metal-rich ejecta of supernovae (SNe) are 
mixed instantaneously and homogeneously into the ISM. 
Models using this approximation, together with theoreti- 
cal nucleosynthesis yields of type la and type II SNe, can 
explain the behaviour of element-to-iron ratios ([el/Fe]) 
of stars as function of metallicity [Fe/H] for many ele- 
ments and for [Fe/H] > —2. This shows that the instan- 
taneous mixing approximation is valid at this stage and - 
since at these metallicities even some of the lowest mass 
core-collapse supernovae (SNe II) have exploded - that 
the stellar yields averaged over the initial mass function 
(IMF) are for most el emen ts accurate within a factor of 
two (see e.g. Samland 1997 ). 

However, observations of very metal-poor stars show 
significant scatter in [el/Fe] ratios at [Fe/H] < —2, im- 
plying that the ISM was not well mixed at this stage 
(Paper I). These local chemical inhomogeneities were 
probably mainly caused by SNe II, since progenitors of 
SN la have much longer lifetimes and are unimportant for 
the chemical enrichment of the ISM until approximately 
[Fe/H] > — 1. At these early stages of galaxy formation, 
the instantaneous mixing approximation is not valid and 
yields depending on the mass of individual SNe II become 
important. Therefore, accurate nucleosynthesis yields as 
a function of progenitor mass are crucial for the under- 
standing of the earliest stages of galaxy formation. 

In Paper I, a stochastic chemical evolution model 
was presented which accounts for local chemical inho- 
mogeneities caused by SNe II with different progenitor 
masses. The model successfully reproduces the scatter in 
[el/Fe] ratios as function of [Fe/H] for some elements like Si 
or Ca, but fails quantitatively in the case of the two most 
abundant a-elements, O and Mg. The scatter in [O/Fc] 
and [Mg/Fe] is much larger than observed and predicts 
stars with [O/Fe] and [Mg/Fe] < -1.0. This result de- 
pends mainly on the employed stellar yields, demonstrat- 
ing that either the oxygen/magnesium or the iron yields 
(or both) as a function of progenitor mass are not well 
determined by existing nucleosynthesis models. 

The solution to this problem is important for the un- 
derstanding of the chemical evolution of our Galaxy. In 
this work, we try to reconcile element abundance observa- 
tions of metal-poor halo stars with the predictions of our 
inhomogeneous chemical evolution model by changing the 
progenitor mass dependence of stellar yields. The forma- 
tion of oxygen and magnesium in hydrostatic burning and 
ejection during a SN event is much better understood than 
the formation and ejection of 56 Ni (which decays to 56 Fe 
and forms the bulk of the ejected iron), since the amount 
of ejected 56 Ni is directly linked to the still not fully under- 



stood explosion m echan ism (c.f. Licbcndo rfcr et al. |2001 



Mezzacappa et al. |200l| ; Rampp & Janka |2000| ). Any at- 
tempt to alter stellar yields should therefore start with 



iron and iron-group elements. We present a method to de- 
rive stellar iron yields as function of progenitor mass from 
the observations of metal-poor halo stars, assuming given 
yields of oxygen and magnesium. 

In Sect. we give a short description of the stochas- 
tic chemical evolution model, followed by a summary of 
observations and basic model results in Sect. 0. The dis- 
cussion of uncertainties in stellar yields and how global 
constraints on stellar iron yields can be gained from ob- 
servations is given in Sect, m Implications for stellar iron 
yields and conclusions are given in Sect. @ and Sect, pi 
respectively. 

2. The chemical evolution model 

Observations of very metal-poor halo stars show a scatter 
in [el/Fe] ratios of order 1 dex. This scatter gradually de- 
creases at higher metallicities until a mean element abun- 
dance is reached which corresponds to the [el/Fe] ratio of 
the stellar yields integrated over the initial mass function. 
Our stochastic chemical evolution model of Paper I follows 
the enrichment history of the halo ISM in a cube with a 
volume of (2.5 kpc) 3 , down to a resolution of (50 pc) 3 . 
Every cell of the grid contains detailed information about 
the enclosed ISM and the mass distribution of stars. For 
the purpose of this paper, the enrichment of the ISM with 
O, Mg, Si, Ca and Fe is computed. 

At every time-step, randomly chosen cells may create 
stars. The likelihood for a cell to form a star is propor- 
tional to the square of the local ISM density. The mass of 
a newly formed star is chosen randomly, with the con- 
dition that the mass distribution of all stars follows a 
Salpeter IMF. The lower and upper mass limits for stars 
are taken to be 0.1 Mq and 50 Mq, respectively. Newly 
born stars inherit the abundance pattern of the ISM out of 
which they formed, carrying therefore information about 
the chemical composition of the ISM at the place and time 
of their birth. 

Stars in a range of 10 — 50 Mq are assumed to explode 
as SNe II (or hypernovae, we will use the term SNe II to in- 
clude hypernovae unless otherwise noted) resulting in an 
enrichment of the neighbouring ISM. Intermediate mass 
stars form planetary nebulae, which return only slightly 
enriched material. Low mass stars do not evolve signifi- 
cantly during the considered time but serve to lock up part 
of the mass, affecting therefore the abundances of elements 
with respect to hydrogen. Stellar yields are taken from 
Thielemann et al. (1996, hereafter TH96) and Nomoto et 
al. (1997). Additionally, since there are no nucleosynthesis 
calculations for 10 Mq progenitors, their yields were set 
to 1/10 of the yields of the 13 Mq model. We then linearly 
interpolated the stellar yields given in these papers, since 
we use a finer mass-grid in our simulation. The interpo- 
lation gives IMF averaged values of [el/Fe] ratios which 
are in good agreement (within 0.1 dex) with the observed 
mean values of metal-poor stars. 

The SN remnant sweeps up the enriched material in 
a spherical, chemically well mixed shell. Since the explo- 
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Fig. 1. [O/Fe] and [Mg/Fe] ratios vs. mctallicity [Fc/H] of metal-poor halo stars (squares and triangles) and model 
stars (dots). Circles depict [O/Fe] and [Mg/Fe] ratios of SN II models of the given progenitor mass. (See text for 
details.) In contrast to observations, model stars with subsolar [O/Fe] and [Mg/Fe] ratios are predicted by the applied 
stellar yields, as visible by the location of the 13 and 15 Mq SNe. 



sion energy of SNe II is believed to depend only sl ightly 
on the mass of its progenitor (Woosley & Weaver 1995, 



hereafter WW95; Thielemann et al. 1996), we assume that 
each SN II sweeps up about 5 x_10_* Mq of gas (Ryan et al. 



1996; Shigeyama & Tsujimoto 1998). Stars which form out 



of material enriched by a single SN inherit its abundance 
ratios and therefore show an element abundance pattern 
which is characteristic for this particular progenitor mass. 
This will lead to a large scatter in the [el/Fe] ratios, as 
long as local inhomogeneities caused by SN events dom- 
inate the halo ISM. As time progresses, supernova rem- 
nants overlap and the abundance pattern in each cell ap- 
proaches the IMF average, leading to a decrease in the 
[el/Fe] scatter at later times. Since the SN remnant for- 
mation is the only dynamical process taken into account, 
this model shows the least possible mixing efficiency for 
the halo ISM. This is just the opposite to chemical evolu- 
tion models which use the instantaneous mixing approxi- 
mation. We continue our calculation up to an average iron 
abundance of [Fe/H] = —1.0. At this metallicity, SN la 
events which are not included in our model start to influ- 
ence the ISM significantly. A more detailed description of 
the model can be found in Paper I. 

We emphasize one important result: Starting with a 
primordial ISM and taking into account local inhomo- 
geneities caused by SNe II, the initial scatter in [el/Fe] 
ratios is determined solely by the adopted nucleosynthesis 
yields. The details of the chemical evolution model only 
determine how fast a chemically homogeneous ISM is 
reached, i.e. how the scatter evolves with time or (equiv- 



alently) iron abundance [Fe/H]. Therefore, the range of 
[el/Fe] ratios of the most metal-poor stars does not de- 
pend on specific model parameters but is already fixed by 
the stellar yields. 



3. Observations and basic model results 



As mentioned in the introduction, existing nucleosynthesis 
models, combined with a chemical evolution model tak- 
ing local inhomogeneities into account, predict [O/Fe] and 
[Mg/Fe] ratios less than solar for some metal-poor stars. 
This is in contrast to observations of metal-poor halo 
stars, as can be seen in Fig. [y. The left hand panel shows 
the [O/Fe] ratio of observed and model stars as func- 
tion of iron abundance [Fe/H] and the right hand panel 
the same for [Mg/Fe], where the model stars are plot- 
ted as small dots. The observational data were collected 
from Magai n (|1989| ), Molaro & Bon ifacio fll99C| ), Molaro 
fc Ca stelli ( |1990D , Peter son et al. fll99CD , Bcssell et a l. 
( |199lD , Ryan et al. (|1991| ), Spiesman fc Wal lerstein ( [1991D , 
Spite & Spite Ql991| ), Norr is et al. fll993| ) Bcv eridge & 
Sneden fl!994|), King (|1994|), Nissen et al. (|1994[), Primas 



et al. (|199J), Sned en et al. (|199j), Fuhrmann et al. (1995), 
Mc William et al. (|1995[ ), Balachandra n fc C arney (1996), 
Ryan et al. (1996), Israelian et al. ( 1998| ), Jchin et al. 
( |1999D , Boesgaa rd et al. ( |1999[ ), Idiart & T hevcn in ( |2000D , 
Carretta et al. (2000) and Israelian et al. ( pOOl ). 

Combining data from various sources is dangerous at 
best, since different investigators use different methods to 
derive element abundances with possibly different and un- 
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known systematic errors. This influences the scatter in 
[el/Fe] ratios, which plays a crucial role in determining 
the chemical (in)homogeneity of the ISM as function of 
[Fe/H]. Unfortunately, there is no investigation with a 
sample of oxygen/magnesium abundances of metal-poor 
halo stars that is large enough for our purpose. Therefore, 
we are forced to combine different data sets, keeping in 
mind that unknown systematic errors can enlarge the in- 
trinsic scatter in element abundances of metal-poor stars. 



Recently, Idiart & Thevenin (200C) and Carretta et al. 
(200C) reanalyzed data previously gathered by other au- 
thors and applied NLTE corrections to O, Mg and Ca 
abundances, which is a first step in reducing the scatter 
introduced by systematic errors. Therefore we divided the 
collected d ata in to two groups, namel y the d ata of Idiart & 
Thevenin ( 2000|) and Carretta et al. ( 2000 ), which is rep- 
resented in Fig. [j] by triangles, and the data of all other 
investigators, represented by squares. If multiple observa- 
tions of a single star exist, abundances are averaged and 
pentagons and diamonds are used for the first and sec- 
ond group, respectively. (Averaging of data points was 
only necessary in a few cases for Mg, Si and Ca abun- 
dances.) Note, that the average random error in element 
abundances is of the order 0.1 dex. 

Also plotted in Fig. [j] as circles are [el/Fe] ratios pre- 
dicted by nucleosynthesis calculations of TH96. The num- 
bers in the circles give the mass of the progenitor star in 
solar masses. In the picture of inhomogeneous chemical 
evolution, a single SN event enriches the primordial ISM 
locally (in our model by mixing with 5 x 10 4 Mq of ISM) 
with its nucleosynthesis products. Depending on the mass 
of the progenitor star, the resulting [O/Fe] and [Mg/Fe] 
ratios in these isolated patches of ISM cover a range of over 
two dex and as long as the ISM is dominated by these local 
inhomogeneities, newly formed stars will show the same 
range in their [el/Fe] ratios. In particular, this means that 
stars with [O/Fe] and [Mg/Fe] as small as —1.0 are in- 
evitably produced by our model. This is in contrast to the 
bulk of observed metal-poor halo stars, which show [O/Fc] 
and [Mg/Fe] ratios in the range between 0.0 and 1.2, and 
is a strong indication that existing nucleosynthesis models 
may correctly account for IMF averaged abundances but 
fail to reproduce stellar yields as function of progenitor 
mass. 



4. Global constraints on stellar Fe yields 

4.1. Uncertainties in O, Mg and Fe yields 

Apart from the shortcomings of nucleosynthesis yields dis- 
cussed in Sect. ||, there seems to be an additional uncer- 
tainty concerning either the stellar yields of O and Mg 
or the derivation of their abundances in metal-poor halo 
stars, as shown in Fig. g: 

The theoretical nucleosynthesis yields of oxygen 
(Yq ( TO )) an d magnesium (Yyig (m)) show a very similar 
dependence on progenitor mass to, i.e. in first order we 
can write Yy[ g (to) « 6.7 ■ 10~ 2 • Yq (to). Thus, for model 
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Fig. 2. [O/Mg] vs. [Mg/H] ratios of metal-poor halo stars. 
Nucleosynthesis models predict a narrow region of possible 
[O/Mg] ratios (hatched) which is not consistent with the 
scatter of observations. Symbols are as in Fig. [y. 

stars [O/Mg] « 0.0 on average, and due to chemical inho- 
mogeneities in the early ISM, model stars scatter in the 
range -0.3 < [O/Mg] < 0.1 (hatched region in Fig. |). In 
contrast to theoretical predictions, observations of metal- 
poor halo stars scatter in the range —0.3 < [O/Mg] < 0.5, 
with a mean of [O/Mg] w 0.15. This result is very impor- 
tant, since it means that either even our understanding of 
nucleosynthesis processes during hydrostatic burning is in- 
complete or that oxygen abundances at very low metallic- 
ities tend to be overestimated (or magnesium abundances 
underestimated) . 

The problem hinted at in Fig. g is also connected to 
the recent finding that the mean [O/Fe] ratio of metal- 
poor halo stars seems to increase with decreasing mctallic- 
ity [Fe/H] , whereas the mean [Mg/Fe] ratio seems to stay 
constant (see e.g. Israelian et al. |199S| , 2001; Boesgaard 
et al. |1999| ; King |2000| but see also Rcbolo et al. [2002p . 
This result can not be explained by changes in the sur- 
face abundances due to rotation, since rotation tends to 
decrease the oxygen abundance in the stellar atmosphere, 
whereas magnesium abundances remain unaffected (Heger 
& Langer |2000| ; Meynet & Maeder |200C| ). However, the 
problem described with Fig. would disappear, if the in- 
crease in [O/Fe] with decreasing metallicity is not real but 
due to some hidden systematic error. Then oxygen abun- 
dances would have to be reduced, resulting in a smaller 
scatter and lower mean in [O/Mg]. 

Regarding nucleosynthesis products, a crude argument 
shows that (at least in the non-rotating case) we should 
not expect a drastic change in the progenitor mass depen- 
dence of O and Mg yields: Oxygen and magnesium are pro- 
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duced mainly during hydrostatic burning in the SN pro- 
genitor and only a small fraction of the ejecta stems from 
explosive neon- and carbon-burning (see e.g. Thielemann 
et al. 



1990 



1996). Magnesium is to first order a prod- 



uct of hydrostatic carbon- and ensuing neon-burning in 
massive stars. The amount of freshly synthesized magne- 
sium depends on the available fuel, i.e. the size of the C-0 
core after hydrostatic He burning, which also determines 
the amount of oxygen that gets expelled in the SN event. 
Thus, O and Mg yields as function of progenitor mass 
should be roughly proportional to each other. A very large 
mass loss during hydrostatic carbon burning could reduce 
the size of the C-0 core and thus decrease the amount of 
synthesized magnesium for a given progenitor mass. This 
would result in a larger scatter of [O/Mg] ratios than in- 
dicated by the hatched region in Fig. y. But the evolu- 
tionary timescale of carbon burning is very short indeed 
(« 5.8 • 10 3 yr for a 25 M star, e.g. Imbriani et al. |200l| ), 
making a significant change in the structure of the C-0 
core unlikely. 

Although the hydrostatic burning phases are thought 
to be well understood, one has to keep in mind that the 
important (effective) 12 C(a, 7) 16 reaction rate is still un- 
certain and that the treatment of rotation and convection 
may also influence the amount of oxygen and magnesium 
produced during hydrostatic burning. Recently, Heger et 
al. (2000) showed that even in the case of slow rotation 
important changes in the internal structure of a massive 
star occur: 

Rotationally induced mixing is important prior to cen- 
tral He ignition. After central He ignition, the timescales 
for rotationally induced mixing become too large com- 
pared to the evolutionary timescales, and the further evo- 
lution of the star is similar to the non-rotating case. 
Nevertheless, He cores of rotating stars are more massive, 
corresponding to He cores of non-rotating stars with about 
25% larger initial mass. Furthermore, for a given mass of 
the He core, the C-0 cores of rotating stars are larger than 
in the non-rotating case. At the end of central He burn- 
ing, fresh He is mixed into the convective core, converting 
carbon into oxygen. Therefore, the carbon abundance in 
the core is decreased, whereas the oxygen abundance is in- 
creased. Unfortunately, no detailed nucleosynthesis yields 
including rotation have been published yet, but since the 
size of the He core is increased in rotating stars, at least 
changes in the yields of a-elements have to be expected. 
(For a review of the changes in the stellar parameters in- 
duced by rotation see Maeder & Meynet 200C .) 

Contrary to oxygen and magnesium which stem from 
hydrostatic burning, iron-peak nuclei are a product of ex- 
plosive silicon-burning. Unfortunately, no self-consistent 
models following the main-sequence evolution, collapse 
and explosion of a massive star exist to date which would 
allow to determine reliably the explosion energy and the 
location of the mass-cut between t he fo rming neutron star 
and t he ejecta (Liebend orfer e t al. 2001 ; Mezzacappa et al. 
2001 ; Rampp & Janka 2000 ). Therefore, nucleosynthesis 
models treat the mass cut usually as one of several free pa- 



rameters and the choice of its value can heavily influence 
the abundance of ejected iron-group nuclei. For this rea- 
son, we feel that oxygen and magnesium yields of nucleo- 
synthesis models are more reliable than iron yields, in spite 
of the uncertainties discussed above. 

To illustrate this point we show a comparison of O 
and Mg yields (Fo (w), Ym s (to), Fig. ||) and of Fe yields 
(Ype (to) , Fig. ^) of nucleosynthesis calculations (neglect- 
ing rotation) from different authors. The models of WW95 
(solar composition "C" models) are marked by filled 
squares, TH96 by filled circles, Nakamura et al. (|200l|, 10 51 



200l| 10 5 ^ 
al. (|2002 



erg models) by open squares and Rauscher et al. (|2002 
"S" models) by open triangles. Upper points in Fig. Hcor 
respond to O yields, lower points to Mg yields. Apart from 
the dip visible in Ym s (m) of the WW95 models, the O and 
Mg yields of the different authors agree remarkably well. 
Chemical evolution calculations by Thomas et al. (1998) 
show that WW95 underestimate the average Mg yield due 
to this dip. The minor differences between the models 
can mostly be attributed to different progenitor models 
prior to core-collapse, the employed 12 C(a, 7) 16 reaction 
rate, the applied convection criterion (e.g. Schwarzschild 
vs. Ledoux) and artificial explosion methods after core- 
collapse (e.g. piston vs. artificially induced shock wave). 
On the other hand, as visible in Fig. [|, Yp (to) of the 
different authors differs by more than an order of magni- 
tude for certain progenitor masses, which is mostly due to 
the arbitrary placement of the "mass-cut" between proto- 
neutron star and ejecta. In order to reconcile the results 
of the inhomogeneous chemical evolution model with ob- 
served [O/Fe] and [Mg/Fe] abundance ratios, it is there- 
fore clearly preferable to artificially adjust Yp e (to) rather 
than Y"o ("m) & n d Fvig (m) . 

For the following discussion, reliable O and Mg yields 
as function of progenitor mass with an estimate of their 
error range are needed. To this end, we calculated best fit 
curves to the To (™) an d Ym s (m) yields of the different 
authors, visible as dashed lines in Fig. g. The low Mg yields 
of the 20 and 22 Mq WW95 models were neglected for 
this purpose. The deviations Ae (to) of the original O and 
Mg yields Y c \ (m) from our best fit yields Y c \ (to) is defined 
as 



Ae (to) 



Y c i (to) - Y c i (to) 
Fei (to) 



(1) 



The error Ae (m) depends on progenitor mass, but 
is generally small. To account for the uncertainty in 
Y e \ (to) introduced by the different nucleosynthesis mod- 
els, we replace in the following the original Y \ (to) by 
(1 ± Ae) • Y e \ (to), where we dropped the dependence of 
Ae (to) on to in the notation. For most progenitor masses, 
Ae < 0.2 for both O and Mg and the maximal deviation is 
in both cases smaller than 0.5. The dotted lines in Fig. g 
show the curves (1 ± Ae) • Y e \ (to) with Ae = 0.2. Since 
Ae is small, the impact of the uncertainties in the stellar 
O and Mg yields is almost negligible for the derivation 
of constraints on Yp c (m). On the other hand, these small 
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Table 1. Best fit O and Mg yields Y c \ (to) proposed in 
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Fig. 3. O and Mg yields from different authors as func- 
tion of progenitor mass. Models are from: WW95 , filled 
squares; TH96, filled circles ; Nak amura et al. ( 2001 ), open 
squares; Rauscher et al. ( |2002 ), open triangles. Upper 
points correspond to O yields, lower points to Mg yields. 
Dashed and dotted lines represent best fit curves to the 
different nucleosynthesis models (see text for details). 
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Fig. 4. Fe yields from different authors as function of pro- 
genitor mass. Symbols are the same as in Fig. [| Contrary 
to O and Mg yields, different authors obtain very different 
yields for a given progenitor mass. This is mostly due to 
the arbitrary placement of the mass cut. 



Sect. 4.1. The first column gives the progenitor mass to 
and the following columns the oxygen and magnesium 
yields (all values in solar masses). 



m 


Yo(m) 


Y Mg (to) 


m 


Fo(to) 


y Mg (to) 


10 


2.2E-02 


1.2E-03 


31 


5.0E+00 


3.3E-01 


11 


8.6E-02 


4.9E-03 


32 


5.4E+00 


3.5E-01 


12 


1.5E-01 


8.5E-03 


33 


5.7E+00 


3.6E-01 


13 


3.1E-01 


2.0E-02 


34 


6.1E+00 


3.8E-01 


11 


4.9E-01 


3.6E-02 


35 


6.5E+00 


4.0E-01 


15 


6.7E-01 


5.3E-02 


36 


6.9E+00 


4.2E-01 


16 


8.7E-01 


6.9E-02 


37 


7.3E+00 


4.4E-01 


17 


1.1E+00 


8.6E-02 


38 


7.7E+00 


4.5E-01 


18 


1.3E+00 


1.0E-01 


39 


8.1E+00 


4.7E-01 


19 


1.5E+00 


1.2E-01 


40 


8.6E+00 


4.9E-01 


20 


1.7E+00 


1.4E-01 


41 


9.0E+00 


5.2E-01 


21 


2.0E+00 


1.5E-01 


42 


9.5E+00 


5.4E-01 


22 


2.2E+00 


1.7E-01 


43 


1.0E+01 


5.6E-01 


23 


2.5E+00 


1.9E-01 


44 


1.0E+01 


5.9E-01 


24 


2.8E+00 


2.0E-01 


45 


1.1E+01 


6.1E-01 


25 


3.0E+00 


2.2E-01 


46 


1.1E+01 


6.4E-01 


26 


3.4E+00 


2.4E-01 


47 


1.2E+01 


6.6E-01 


27 


3.7E+00 


2.6E-01 


48 


1.2E+01 


6.8E-01 


28 


4.0E+00 


2.7E-01 


49 


1.3E+01 


7.1E-01 


29 


4.3E+00 


2.9E-01 


50 


1.3E+01 


7.3E-01 


30 


4.6E+00 


3.1E-01 









uncertainties may (almost) be able to explain the discrep- 
ancy in the scatter of [O/Mg] ratios between observations 
and model stars visible in Fig. 0. Allowing for a mean devi- 
ation of Ae = 0.2, the maximal scatter in [O/Mg] over all 
progenitor masses may extend to —0.25 < [O/Mg] < 0.3, 
which is very close to the one observed. Future nucleo- 
synthesis calculations have to show whether this interpre- 
tation is correct or not. 

For the remaining discussion, we adopt the term 
(1 ± Ae) • Y e \ (to) with Ae = 0.2 for the stellar oxygen 
and magnesium yields as a function of progenitor mass, 
assuming that they reproduce the true production in mas- 
sive stars well enough. The values adopted for the best fit 
yields Y e \ (to) are given in Table [y. 

4.2. The influence of Z and SNe from Population III 
stars 

Apart from the uncertainties in the O and Mg yields dis 
cussed in Sect. 



4.1 



, nucleosynthesis yields may also depend 
on the metallicity of the progenitor. Unfortunately, the 
question how important metallicity effects are is far from 
solved. Nucleosynthesis calculations in general neglect ef- 
fects of mass loss due to stellar winds. WW95 present 
nucleosynthesis results for a grid of metallicities from 
metal-free to solar and predict a decrease in the ejected 
O and Mg mass with decreasing metallicity. However, the 
O yields presented lie all in the range covered by the best 
fit yields (1 ± Ae) • Y c \ (to) with Ae = 0.2 adopted for 
this paper. This is not true in the case of Mg. But since 
the "dip" in Y^ g (to) visible in Fig. [| gets more and more 
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pronounced with lower Z and since it is known (Thomas 
et al., 1998), that WW95 underestimate the average pro- 



duction of Mg, we feel that the metallicity dependence of 
Mg yields is not established well enough to include this 
feature into our analysis. 

Co ntrary to the results obtained by WW95, Maeder 
( |l992|) showed that stellar O yields decrease with increas- 
ing metallicity due to strong mass loss in stellar winds. 
(No detailed results were given in the case of magne- 
sium.) Stars more massive than 25Mq with solar metal- 
licity (Z=0.02), eject large amounts of He and C in stel- 
lar winds (prior to the conversion into oxygen) which re- 
sults in dramatically reduced O yields. Metal-poor stars 
(Z<0.001) do not undergo an extended mass-loss phase 
and their O yields are comparable to the ones given by 
WW95 and TH96. Since Z<0.001 roughly corresponds to 
[Fe/H] < —1.5 and we are mainly concerned with metal- 
poor halo stars in this metallicity range, we can neglect 
changes in O yields due to metallicity. 

Recently, Heger & Woosley ( 2000| ) published nucleo- 
synthesis calculations of pair-instability SNe from very 
massive, metal-free (Population III) stars in the mass 
range from 140Mq to 260Mq. For Population III stars in 
the mass range 25 — 140Mq and stars more massive than 
260Mq, black hole formation without ejection of nuc leo- 
synthesis yields seems likely (Heger & Woosley 2000| ). In 
order to investigate the influence of those massive metal- 
free stars on the enrichment of the ISM and especially 
their impact on the distribution of model stars in [O/Mg] 
vs. [Mg/H] plots (c.f. Sect. 4.1), we carried out several in- 
homogeneous chemical evolution calculations with vary- 
ing SF efficiencies and IMF shapes for the Population III 
stars. The detailed results will not be shown here, but 
some basic conclusions are discussed in the following: 

The theoretical scatter in [O/Mg] predicted by the 
massive Population III stars lies in the range —0.3 < 
[O/Mg] < 0.3. This could help to explain the scatter 
in [O/Mg] observed in metal-poor stars, if we take ob- 
servational errors of the order of 0.1 dex into account. 
Models with a high SF efficiency for Population III stars 
indeed show some stars with high [O/Mg] ratios. However, 
61% of the metal-poor stars ([Fe/H] < —1.0) with ob- 
served O and Mg abundances show [O/Mg] > 0.1 (see 
Fig. g), whereas less than 1% of the model stars have 
[O/Mg] ratios in this range (the exact number depends 
on the shape of the IMF). Clearly, the observations of 
metal-poor stars can not be explained as a consequence 
of such massive, metal-free SNe. Furthermore, the distri- 
bution of model-stars in [O/Fe] and [Mg/Fe] can not be 
reconciled with observations of metal-poor stars. If the SF 
efficiency of Population III stars is small, these discrepan- 
cies in [O/Fe] and [Mg/Fe] disappear, but the number of 
model stars with [O/Mg] > 0.1 decreases even further. We 
therefore conclude that - at least for the purpose of this 
paper - the (possible) influence of SNe from very massive 
Population III stars can safely be neglected. 



4.3. Putting constraints on Fe yields with the help of 
observations 

In order to reproduce the scatter of observed [O/Fe] and 
[Mg/Fe] ratios of metal-poor halo stars while keeping the 
oxygen and magnesium yields fixed, we have to adjust the 
stellar iron yields Yp c (to) as function of progenitor mass 
m. Since it is not known from theory what functional form 
Yp G (to) follows (increasing with m, declining or a more 
complex behaviour), we have the freedom to make some 
ad hoc assumptions. Nevertheless, some important con- 
straints on Yp e (to) can be drawn from the scatter, range 
and mean of observed [O/Fe] and [Mg/Fe] abundances, as 
visible from Fig. Qt 

1. IMF averaged stellar yields (integrated over a complete 
generation of stars) should reproduce the mean oxygen 
and magnesium abundances of metal-poor halo stars, 
i.e. [O/Fe] w 0.4 and [Mg/Fe] w 0.4. 

2. Stellar yields have to reproduce the range and scatter 
of [el/Fe] ratios observed. Using oxygen and magne- 
sium as reference this requires: 



0.0 < [O/Fe] < 1.2, 
0.1 < [Mg/Fe] < 1.2. 



(2) 
(3) 



(Note, that the error in abundance determinations is 
of order 0.1 dex.) 

There exist a few Type II and Type Ib/c SN ob- 
servations (1987A, 1993J, 19941, 1997D, 1997ef and 
1998bw) where the ejected 56 Ni mass and the mass 
of the progenitor was derived by analyzing and mod 
elling the light-curve ( e.g. S untzeff & Bouchet 
Shigeyama & Nomoto 1990; Shigeyama et al 



1990 



1994 



Iwam oto et al. |l994| , h998 [ [2000^ Kozma & Fra nsson 
1998| ; Turatto et al. L99§ Chugai & Utrobin [2000 ; 
Sollerman et al. 200C). These observations give im- 
portant constraints on Yp e (to) since they constrain 
the stellar yields for some progenitor masses, although 
they are not unambiguous (see Sect. 4.3.3| ). 
Since observations of metal-poor halo stars show no 
clear trends in [Mg/Fe] with decreasing [Fe/H] we re- 
quire that modified iron yields likewise do not intro- 
duce any skewness in the distribution of model stars. 
In the case of oxygen, it is not clear yet whether the 
apparent slope in [O/Fe] in recent abundance studies 
is real or due to some systematic errors (see Fig. Hand 
Sect. O). 



It is clear that it is not possible to predict lpc (to) 
unambiguously, since the information drawn from obser- 
vations is afflicted by errors. We therefore do not attempt 
to find a solution which reproduces the observations per- 
fectly, but try to extract the global properties of Yp e (to) 
needed to explain the behaviour of observed [el/Fe] ratios 
in metal-poor halo stars. 
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4.3.1. IMF averaged iron yields 

Since the yields of TH96 were calibrated so that the first 
constraint is fulfilled, we require that the average iron 
yield of SNe II stays constant when we change the pro- 
genitor mass dependence of Y" Fc (to). Assuming a Salpeter 
IMF ranging from 0.1 Mq to 50 Mq and assuming that 
all stars more massive than 10 Mq turn into core -colla pse 
SNe (or hypernovae, see e.g. Nakamura et al. 2001 ), a 
SN II produces on average (Yo) ~ 1.9 Mq of oxygen, 
(F Mg ) ~ 0.12 Mq of magnesium and (r Fc ) w 0.095 M Q 
of iron. Leaving the average oxygen and magnesium yields 
unchanged and modifying the stellar iron yields, we there- 
fore always have to require that on average « 0.095 Mq of 
iron are ejected per SN. Thus, Y Fc (to) has to satisfy the 
following condition: 



.Ml 



(Y Fe ) - 



Yp e (to) to" 



-2.35 



dm 



u 



50 
10 '' 



2 - 35 dm 



0.095 M Q . 



(4) 



Note that the average [el/Fe] ratio of the model stars 
depends on the lower and upper mass limits of stars that 
turn into SNe II and their yields. If we raise the lower 
mass limit, the average oxygen yield of a SN will increase 
since many stars with a low oxygen yield no longer con- 
tribute to the enrichment of the ISM. The same is true for 
magnesium and iron and the combination of the new av- 
eraged yields may lead to slightly changed average [el/Fe] 
ratios. Since there are only a few SNe with large progeni- 
tor masses, changing the upper mass limit will have a very 
small influence on the average [el/Fe] ratios. For the re- 
maining discussion, we will keep the lower and upper mass 
limits of SNe II fixed at 10 Mq and 50 Mq, respectively. 

4.3.2. Range and scatter of observations 

The second constraint can be used to calculate the range 
that Yp e {m) has to cover. In our picture of inhomoge- 
neous chemical evolution, we assume that the first SNe 
locally enrich the primordial ISM. Stars forming out of 
this enriched material therefore inherit the [el/Fe] ratios 
produced by these SNe which is determined in turn by 
the stellar yields Y c \ (to). (For the time being, we neglect 
the additional uncertainty hidden in the factor (1 ± Ae).) 
Thus, for the first few generations of stars formed at the 
time the ISM is dominated by local chemical inhomo- 
geneities, the following identity holds (with the exception 
of H and He where also the abundances in the primordial 
ISM have to be taken into account): 

[el/Fe] = log — — /Ar — = log 



log 



F el (to) /Y Fe (m) 



where iV cl q (N e \^) is the number density of a given ele- 
ment el in the solar (stellar) atmosphere and M ol q (M e i ]Vt ) 
the corresponding mass fraction. (Solar abundances were 



taken from Anders & Grevesse 198£). Now, let Y e \ (to) be 
either the stellar yields of oxygen or of magnesium and 
a, /3 the minimal and maximal [cl/Fc] ratios derived from 
observations. Then the constraint gives: 

a < [el/Fe] < ^=^ 



Y_^)JY^ < 
M clB /M FcB -" 



Y cl (m) • 10-' 3 • -jp@. < y Fc ( m ) 

M el,0 



< Fel (TO) ' 10" 



M Fe,Q 



M, 



1,0 



Since the yields of oxygen and magnesium as function of 
progenitor mass are assumed to be known, we now have 
two sets of inequalities for the stellar iron yields. The first 
is derived from the minimal and maximal [O/Fe] ratios 
(Eq.|): 

8.37 • 10~ 3 • Y (to) < y Fc (m) < 1.33 • 10" 1 • Y (to) , (5) 

and the second from the minimal and maximal [Mg/Fe] 
ratios (Eq. ||): 

1.23 • 10" 1 • Y Mg (to) < Y Fc (to) < 2.46 ■ F Mg (to) , (6) 

where the uncertainty in the O and Mg yields given by the 
factor (1 ± Ae) was neglected. Thus, for any given progen- 
itor mass, y Fe (to) is only determined within a factor of 20 
- 25 and further constraints are needed to derive a reliable 
iron yield. 

Fig. g shows the stell ar iro n yields Y Fc (to) (soli d line ) 
from Thielemann et al. ( 1996] ) and Nomoto et al. ( 1997 ), 
binned with a bin size of 1 Mq. To reproduce the range 
and scatter of [O/Fe] and [Mg/Fe] ratios observed in 
metal-poor halo stars, Y Fc (m) should remain in the re- 
gion enclosed by the boundaries given by Eqs. (^|) and (|6|), 
which are shown as dashed (representing [O/Fe] = 0.0), 
long dashed (representing [O/Fe] = 1.2), dotted (repre- 
senting [Mg/Fe] = —0.1) and dash-dotted (representing 
[Mg/Fe] = 1.2) lines. Note that the lower lines in Fig. | 
represent the upper boundaries derived from metal-poor 
halo stars and vice versa. Evidently, the iron yields of SNe 
with progenitor masses in the range 10 — 15 Mq are out- 
side the boundaries given by metal-poor halo stars, leading 
to the stars in our model with much too low [O/Fe] and 
[Mg/Fe] abundances (Paper I, c.f. also Fig. pi). Therefore, 
we can already conclude that the iron yields of these SNe 
have to be reduced to be consistent with observations. 
Consequently, the iron yields of some higher-mass SNe 
have to be increased to keep the IMF averaged [el/Fe] 
ratios constant (Eq. (Q)). This can easily be achieved by 
assuming a higher explosion energy than the "canonical" 
10 51 erg of standard SN models for the more massive stars 
(M > 30 Mq), as was shown recently by Nakamura et al. 
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[Mg/Fe] =1.20 - 

[Mg/Fe] = -0-10 
[O/Fe] = 1.20 
[O/Fe] = 0.00 



suits agree remarkably well. The progenitor mass was es- 
timated to be 20 ± 2M , while O.O7O±^M of 56 Ni 



20 30 40 

Progenitor Main Sequence Moss [Msol] 



50 



Fig. 5. Ype (tti) as function of progenitor mass and bound- 
aries constraining the range stellar iron yields have to sat- 
isfy to reproduce the scatter in [O/Fe] and [Mg/Fe] of 
metal-poor halo stars. A ccordi ng to nucleosynthesi s mod - 
els of Thielemann et al. ( |l996[ ) and Nomoto et al. ( |l997j ), 
SNe in the range 10 — 15 Mq clearly eject too much iron 
to be consistent with observations. Also shown are obser- 
vations of core-collapse SNe with known progenitor mass 
and ejected 56 Ni mass. 



(2001). The reader should note that Thielemann et al. 



(1996) and Nomoto et al. (1997) calculated models only 
for the 13, 15, 18, 20, 25, 40 and 70 M progenitors. For 
the 10 M progenitor we assumed an ad hoc iron yield of 
one tenth of the yield of a 13 M star and interpolated 
the intermediate data points. The details of the interpo- 
lation and especially the extrapolation down to the 10 
M star influence the mean [el/Fe] ratio of the ISM at 
late stages, when it can be considered chemically homoge- 
neous. However, this does not change the conclusion that 
the 13 and 15 M models of TH96 produce too much 
iron (if we assume the oxygen and magnesium yields to 
be correct), as is evident from Figs. |] and ||. 



4.3.3. 56 Ni yields from observed core-collapse SNe 

There are six core-collapse SNe with known progenitor 
mass and ejected 56 Ni mass (which is the main source 
of 56 Fe in SNe II explosions, by the decay 56 Ni ->■ 56 Co 
-> 56 Fe), namely 1987A, 1993J, 19941, 1997D, 1997ef and 
1998bw, that are shown in Fig. |. Of these, SN 1987A is 
the m ost extensively studied (see e .g. Suntzcff & B ouchet 
1990J; Shigeyama & N omoto |l990| ; Bouchet et al. |l991a , 



1991b: Suntzcff et al. 1992; Bouchet fc Danziger 1993 



Kozma & Fransson 1998; Fryer et al. 1999) and the re 



were ejected during the SN event. 

SN 1993J had a progenitor in the mass range between 
12 to 15 M and ejecte d approximately 0.0 8 Mp- of 56 Ni 
(Shigeyama et al. 1994; Houck & Fransson 1996), which 
is very similar to SN 19941 with a 13 to 15 M progeni- 
tor and 0.075 M of ejected 56 Ni (Iwamoto ct al. |1994[ ). 
Although the amount of 56 Ni ejected by those SNe lies 
at the upper limit allowed under our assumptions, Fig. a 
shows that these values are still consistent with the con- 
straints given by Eqs. g and g. 

Also consistent with our constraints is SN 1997ef with 
a progenitor mass of 30 — 35 M^ and a 56 Ni mass of 



0.15 ± 0.03 M (Iwamoto et al. gOOGJ). The correspond- 
ing iron yield of SN 1997ef is higher than predicted by 
the nucleosynthesis models of TH96. This is exactly the 
behaviour needed to adjust Yp e (m) according to our con- 
straints. SN 1997ef does not seem to be an ordinary core- 
collapse supernova. The model with the best fit to the 
lightcurve has an explosion energy which is about eight 
times higher than the typical 10 51 erg of standard SN mod- 
els. Such hyperenergetic Type Ib/c SNe are also termed 
hypernovae and probably indicate a change in the explo- 
sion mechanism around 25 — 30 M which could result in 
a discontinuity in the iron yields in this mass range. 

In the case of SN 1997D, the situation is not clear. 
Turatto et al. (1998) propose two possible mass ranges 
for its progenitor: They favour a 26 M star (although 
the progenitor mass can vary from 25 — 40 M ) over a 
possible 8—10 M progenitor. A recent investigation by 
Chugai & Utrobin ( |2000 ) implies a progenitor mass in the 
range 8 — 12 Mq. Both groups deduce an extremely small 
amount of newly synthesized 56 Ni of only « 0.002 M and 
an unusual low explosion energy of only a few times 10 
erg. Since the situation about the progenitor mass remains 
unclear, both possible mass ranges are shown in Fig. ^[ 
On the basis of the small amount of synthesized oxygen of 
only 0.02- 0.07 M (Chugai & Utrobin |200(i| ) , we strongly 
favour the low-mass progenitor hypothesis, since accord- 
ing to nucleosynthesis calculations by TH96 and WW95 
a high-mass progenitor would produce a large amount of 
oxygen (« 3 M for a 25 M progenitor). Moreover, in the 
latter case the observed 56 Ni abundance lies completely 
outside the boundaries derived in Sect. 4.3.2, as can be 
seen in Fig. [| 

SN 1998bw seems to be another hypernova with a ki- 
netic energy of (2 — 5) x 10 52 erg and may be physically 
connected to the underluminous 7-ray burst GRB980425 
(Galama et al. 1998| ; Iwamoto et al. 1998; Iwamoto 1999a, 
|l999b ) . The hypernova model assumes a progenitor mass 



of about 40 M , ejecting w 0.7 M of 5b Ni. Recently. 



Sollerman et al. ( [2000 ) observed SN 1998bw at late phases 
and made detailed models of its light curve and spectra. 
They propose two possible scenarios for this hypernova: 
one with a progenitor mass of 40 M and ejected nickel 
mass of 0.5 M and the other with a 25 M progenitor 
and 0.9 M of nickel. Note that the amount of nickel pre- 
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Fig. 6. Iron yields Yp c (to) respecting the constraints given 
by observations of metal-poor halo stars. Ype (to) starts 
at a very low value and increases continuously. A linear 
increase is not possible since the mean [el/Fe] ratios have 
to be conserved. 



sumably synthesized by this 25 Mq SN is about 10 times 
larger than predicted by SN models that use the "canoni- 
cal" kinetic explosion energy of 10 51 erg. Nevertheless , it is 
still consistent with the constraints derived in Sect. 4.3.2 



and with recent calculations of explosive nucleosynthesis 
in hypernovae by Nakamura et al. (2001). 



4.3.4. Slopes in [el/Fe] vs. [Fe/H] distributions 



Using only the constraints discussed in Sects. 4.3.1 and 



4.3.2 still allows for a wide variety of possible iron yields 
Yfc (to). This is demonstrated by Figs. gand||, where two 
simple ad hoc iron yield functions are shown. The distri- 
butions of model stars resulting from these iron yields are 
plotted in Figs. |7J and ||. 

In Fig. ^ the iron yield Yp c (to) starts at the [O/Fe] 
= 1.2 boundary, increases continuously with increasing 
progenitor mass and ends at the [O/Fe] — 0.0 boundary. 
Consequently, low-mass SNe create a high [O/Fe] ratio in 
their surrounding primordial ISM, whereas it is close to so- 
lar in the neighbourhood of high-mass SNe. The resulting 
[O/Fe] distribution of model stars can be seen in Fig. M. 
The distribution shows a clear trend from high to low 
[O/Fe] ratios with increasing [Fe/H]. A simple least-square 
fit to our data yields [O/Fe] = -0.21 x [Fc/H] +0.01. This 
i s in s urprisingly good agreement with the result of King 
(|2000|) , who finds the relation [O/Fe] = -0.18 x [Fe/H] + 
0.02 after considering the effects of NLTE corrections to 
oxygen abundance determinations from UV OH-lincs. 



Fig. 7. Distribution of [O/Fe] ratios vs. metallicity [Fe/H] 
of model stars resulting from the iron yields shown in 
Fig. @. The slope in the distribution of the model stars is 
consistent with observations of oxygen abundances. A sim- 
ilar slope is introduced in the [Mg/Fe] distribution which 
can not be reconciled with observed magnesium abun- 
dances (See text for details, symbols are as in Fig. |l|). 
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Fig. 8. Iron yields Ifc (to) respecting the constraints given 
by observations of metal-poor halo stars. Yp e (to) starts at 
a somewhat higher value than in Fig. |6|, reaches a maxi- 
mum at about 30 Mq and decreases again. 
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Fig. 9. Distribution of [O/Fe] ratios vs. metallicity [Fe/H] 
of model stars resulting from the iron yields shown in 
Fig. g. Clearly, the rising slope in the distribution of model 
stars is not consistent with observations of metal-poor halo 
stars. Symbols are as in Fig. [I]. 

The reason for the slope in our model is given by the 
distribution of [O/Fe] and [Fe/H] ratios induced in the 
metal-poor ISM by core-collapse SNe with different pro- 
genitor masses, as indicated by the position of the cir- 
cles in Fig. 0. All the low mass SNe with progenitors up 
to 15 Mq induce high [O/Fe] and low [Fe/H] ratios in 
the neighbouring primordial ISM, whereas high-mass SNe 
produce low [O/Fe] and high [Fe/H] ratios (recall that 
a constant mixing mass of 5 x 10 4 Mq per SN event is 
assumed, c.f. Sect. |2|). SNe with intermediate masses in- 
duce [O/Fe] ratios that lie approximately on a straight line 
connecting the two extrema. The distribution of model 
stars of the first few stellar generations follows this line 
closely. As the mixing and chemical enrichment of the 
halo ISM proceeds, the distribution then converges to the 
IMF averaged [O/Fe] ratio. Although the inhomogeneous 
enrichment is responsible for the slope which is in good 
agreement with observations by e.g. Israelian et al. ( 199a ), 
Boesgaard et al. ( p.999[ ), and King ( |2000| ), it fails to re- 
produce the scatter seen in observed oxygen abundances. 
Model stars with [O/Fe] w 1.2 exist only at [Fe/H] < -3.5 
and not at [Fe/H] rs —2.5, where several are observed. 
Furthermore, a similar slope is introduced in the [Mg/Fe] 
distribution ([Mg/Fe] = -0.26 x [Fe/H] - 0.07), where 
none is seen in observations and several model stars show 
[Mg/Fe] ratios as large as ~ 1.5. Therefore, th is Yp c (to) 
has to be rejected (see however Rebolo et al. 2002 con- 
cerning [Mg/Fe]). 

The situation displayed in Figs. || and [| is even worse. 
Here, the iron yield starts at the [O/Fe] = 0.0 bound- 



ary, increases with progenitor mass, reaches its maximum 
around 30 Mq, decreases again and ends at the [O/Fc] 
= 1.2 boundary. The resulting distribution in [O/Fc] 
shows a rising slope, which is clearly in contradiction with 
observations. The slope is a consequence of the fact that 
the iron yields Yp e (jn) of SNe in the range 10 — 15 Mq 
stay very close to the boundary that represents the ratio 
[O/Fe] = 0.0. SNe in this mass range form the bulk of 
SN II events and it is therefore not surprising that their 
large number introduces such a slope in the distribution of 
model stars. Thus, this Yp c (to) also has to be discarded. 
These simple examples show that Ir, (m) should not 
run parallel to the boundaries over a large progenitor mass 
interval, otherwise an unrealistic slope is introduced in 
the [el/Fe] distribution of model stars. They demonstrate 
further, that the information drawn from metal-poor halo 
stars alone is not sufficient to derive reliable iron yields, 
and that information from SN II events and the shape of 
the [el/Fe] distribution as function of [Fe/H] (i.e. how fast 
the scatter decreases and whether slopes are present or 
not) has to be included in our analysis. 

5. Implications for stellar Fe yields 



In Sect. (4.3.4) we have shown that the iron yields of 
lower-mass SNe (in the range 10 — 20 Mq) are crucial to 
the distribution of model stars in [el/Fe] vs. [Fe/H] plots, 
since progenitors in this mass range compose the bulk (ap- 
proximately 69%) of SN II events. Thus, the iron yields of 
lower-mass SNe should not introduce a slope in the [el/Fe] 
distribution but should cover the entire range of observed 
[O/Fe] and [Mg/Fe] ratios in order to reproduce the ob- 
servations. To accomplish this, Yp c (to) should not lie too 
close to the boundaries given in Eqs. (§) and (g) in this 
mass range but should start at the lower boundary ([O/Fe] 
= 1.2) and reach the upper boundary ([O/Fc] = 0.0) for 
some progenitor in the mass range 10 — 20 Mq. If the ob- 
served 56 Ni production of SN 1993J, 19941 and 1987A are 
also taken into account, the observational constraints are 
stringent enough to fix the iron yields of the low mass SNe 
apart from small variations: Yp e (m) starts at the lower 
boundary, increases steeply in the range 10 — 15 Mq to 
the values given by SN 1993J and 19941 and remains al- 
most constant in the range 15 — 20 Mq (to account for SN 
1987A). For the remaining discussion we therefore assume 
the Fe yields in this range to be w 1.5 • 10~ 4 Mq for a 
10 Mq progenitor, ps 5.5 • 10~ 2 Mq for a 15 Mq progen- 
itor and w 7.0 • 10~ 2 M Q for a 20 M progenitor. (The 
detailed yields resulting from our analysis are listed in 
Table 0). We now take a look at possible iron yields of 
higher mass SNe corresponding to the different models of 
the progenitor masses of SN 1997D and 1998bw. There 
are four possible combinations of the progenitor masses of 
those two SNe. 

SI: The first case (model SI, shown in Fig. |10|) gives 
the best fit to abundance observations of metal-poor halo 
stars. Here, we preferred the lower mass progenitor mod- 
els of SN1997D and SN 1998bw over the higher mass 



12 



D. Argast et al.: Implications of O and Mg abundances for stellar iron yields 




[Mg/Fe] = 1.20 
[Mg/Fe] = -0-10 
[O/Fe] = 1.20 
[O/Fe] = 0.00 




[Mg/Fe] = 1.20 
[Mg/Fe] = -0.10 
[O/Fe] = 1.20 
[O/Fe] = 0.00 



20 30 40 

Progenitor Moin Sequence Moss [Msol] 



20 30 40 

Progenitor Main Sequence Mass [Msol] 



Fig. 10. Model SI: Iron yields Yp B (m) respecting the con- Fig. 12. Model H2: Iron yields Yp e (m) respecting the con- 
straints deduced from metal-poor halo stars and SN obser- straints deduced from metal-poor halo stars and SN obser- 
vations. The 10 M model of SN 1997D and 25 M model vations. The 26 M model of SN 1997D and 25 M model 
of SN 1998bw arc assumed to be correct. of SN 1998bw are assumed to be correct. 



models. The curve is characterized by a peak of 0.59 M served in metal-poor halo stars. Obvi ously, l p c (m ) ful- 

of iron at 25 M and a slow decline of the yields down fils th e constraints discussed in Sects. ( 4.3.1 ), ( 4.3.2 ) and 

to 9.5 ■ 1O _2 M for the 50 M progenitor. The yields ( 4.3.3| ). Since no slope is visible in the resulting distribu- 

have to decline again to meet the mean abundances ob- tion of [O/Fe], [Mg/Fe], [Si/Fe] and [Ca/Fe] ratios (shown 
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Fig. 11. Model HI: Iron yields Yp e (m) respecting the con- Fig. 13. Model S2: Iron yields Yp e (m) respecting the con- 
straints deduced from metal-poor halo stars and SN obser- straints deduced from metal-poor halo stars and SN obser- 
vations. The 26 M model of SN 1997D and 40 M model vations. The 10 M model of SN 1997D and 40 M model 
of SN 1998bw are assumed to be correct. of SN 1998bw are assumed to be correct. 
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Table 2. Iron yields Yp c (to) proposed in Sect. 0. The first 
column gives the progenitor mass to in solar masses. The 
following columns give the iron mass (in solar masses) 
synthesized in the corresponding SN event according to 
nucleosynthesis calculations of Thielemann et al. (1996) 
and Nomoto et al. ( [1997D - denoted by TN - and the 
models SI, S2, HI and H2. 



TN 



SI 



S2 



HI 



H2 



10 


1.6E-02 


1.5E-04 


1.6E-04 


1.5E-04 


1.5E-04 


11 


8.0E-02 


1.9E-03 


1.9E-03 


1.8E-03 


1.8E-03 


12 


1.3E-01 


6.8E-03 


6.8E-03 


6.7E-03 


6.8E-03 


13 


1.6E-01 


2.0E-02 


2.0E-02 


1.9E-02 


1.9E-02 


14 


1.6E-01 


5.2E-02 


5.2E-02 


5.1E-02 


5.1E-02 


15 


1.4E-01 


5.5E-02 


5.5E-02 


5.6E-02 


5.6E-02 


16 


1.2E-01 


6.0E-02 


6.0E-02 


6.0E-02 


6.0E-02 


17 


1.1E-01 


6.3E-02 


6.3E-02 


6.3E-02 


6.3E-02 


18 


1.0E-01 


6.7E-02 


6.7E-02 


6.7E-02 


6.7E-02 


19 


9.0E-02 


7.0E-02 


6.9E-02 


6.9E-02 


7.1E-02 


20 


8.0E-02 


7.1E-02 


7.0E-02 


6.9E-02 


7.1E-02 


21 


7.4E-02 


1.1E-01 


7.0E-02 


7.1E-02 


1.1E-01 


22 


7.0E-02 


1.8E-01 


7.0E-02 


7.3E-02 


1.8E-01 


23 


6.5E-02 


2.7E-01 


7.0E-02 


7.6E-02 


2.7E-01 


24 


6.1E-02 


4.2E-01 


7.0E-02 


7.8E-02 


4.2E-01 


25 


5.8E-02 


5.9E-01 


7.0E-02 


8.1E-02 


5.9E-01 


26 


5.9E-02 


5.1E-01 


7.0E-02 


1.9E-03 


1.9E-03 


27 


6.1E-02 


4.3E-01 


7.0E-02 


9.6E-03 


6.7E-03 


28 


6.3E-02 


3.7E-01 


7.0E-02 


2.3E-02 


1.6E-02 


29 


6.5E-02 


3.2E-01 


7.0E-02 


4.9E-02 


3.4E-02 


30 


6.7E-02 


2.7E-01 


8.2E-02 


8.2E-02 


5.8E-02 


31 


6.9E-02 


2.2E-01 


9.8E-02 


1.3E-01 


8.9E-02 


32 


7.1E-02 


2.0E-01 


1.2E-01 


1.8E-01 


1.3E-01 


33 


7.3E-02 


1.8E-01 


1.4E-01 


2.4E-01 


1.7E-01 


34 


7.5E-02 


1.6E-01 


1.7E-01 


3.0E-01 


2.1E-01 


35 


7.7E-02 


1.5E-01 


2.0E-01 


3.5E-01 


2.4E-01 


36 


7.9E-02 


1.4E-01 


2.3E-01 


3.9E-01 


2.7E-01 


37 


8.1E-02 


1.3E-01 


2.8E-01 


4.4E-01 


3.1E-01 


38 


8.3E-02 


1.3E-01 


3.3E-01 


4.8E-01 


3.3E-01 


39 


8.5E-02 


1.2E-01 


3.9E-01 


5.3E-01 


3.6E-01 


40 


8.6E-02 


1.2E-01 


4.6E-01 


5.8E-01 


3.9E-01 


41 


8.7E-02 


1.1E-01 


5.3E-01 


6.4E-01 


4.3E-01 


42 


8.7E-02 


1.1E-01 


6.3E-01 


6.9E-01 


4.5E-01 


43 


8.7E-02 


1.0E-01 


7.6E-01 


7.7E-01 


4.9E-01 


44 


8.7E-02 


1.0E-01 


8.8E-01 


8.6E-01 


5.3E-01 


45 


8.7E-02 


9.8E-02 


1.0E-00 


9.3E-01 


5.8E-01 


46 


8.7E-02 


9.6E-02 


1.2E-00 


1.0E-00 


6.1E-01 


47 


8.7E-02 


9.6E-02 


1.4E-00 


1.1E-00 


6.5E-01 


48 


8.7E-02 


9.5E-02 


1.6E-00 


1.2E-00 


7.0E-01 


49 


8.6E-02 


9.5E-02 


1.8E-00 


1.2E-00 


7.4E-01 


50 


8.6E-02 


9.5E-02 


2.1E-00 


1.3E-00 


7.9E-01 



in Fig. Oh, the constraint described in Sect. ( |4.3.4 ) is 
also respected. The distribution of model stars in [O/Fe], 
[Mg/Fe] and [Si/Fe] is in good agreement with the distri- 
bution of observed stars, whereas a few stars with too low 
[Ca/Fe] ratios are predicted. However, this may be due 
to the fact that Ca stems from explosive O and Si burn- 
ing and therefore depend on the structure of the progeni- 
tor model and the (assumed) explosion energy (Paper I). 
Note, that the mean [Mg/Fe] and [Ca/Fe] ratios in the 



[el/Fe] plots are slightly shifted compared to the mean of 
observations. This problem also occurs when the original 
yields of TH96 are used (Paper I) and will persist for ev- 
ery Ipe (to) we present, since we did not change the mean 
iron yield of 0.095 Mq (Eq. (0)). Especially noteworthy is 
the good agreement in [Si/Fe], since we did not include 
Si in the derivation of the constraints discussed above. 
Moreover, the hypothetical iron yields in the other mod- 
els below all result in [Si/Fe] distributions that do not fit 
the observations as well as model SI. Therefore, we feel 
that model SI gives the best ht to element abundances in 
metal-poor halo stars. 

HI: Fig. [by shows the iron yields under the assumption 
that the higher mass models of SN 1997D and SN 1998bw 
are correct (model HI). Here Yp e (to) stays almost con- 
stant up to 25 Mq, followed by a sudden plunge of the 
yields down to 1.9 • 1O _3 M to account for SN 1997D 
and then a continuous rise to 0.79 Mq of synthesized iron 
for the 50 Mq progenitor that is necessary to account for 
the IMF averaged yield. This sudden decrease of the iron 
yields could indicate a change in the explosion mecha- 
nism from supernovae with "canonical" kinetic explosion 
energies of 10 51 erg to hypernovae with 10 — 100 times 
higher explosion energies. However, as visible in Fig. nil 
the very low iron yield of the 26 Mq progenitor violates 
the [O/Fc] = 1.2 and [Mg/Fe] = 1.2 boundaries derived 
from observations, so we would expect model stars with 
much too high [O/Fe] and [Mg/Fe] ratios. This is indeed 
the case, as can be seen in the corresponding [el/Fe] dis- 
tribution (Fig. nS). A closer examination of the [el/Fe] 
distributions, on the other hand, reveals that these model 
stars are mainly present at very low metallicities ([Fe/H] 

< —2.5). This makes it difficult to decide, whether model 
HI with its dip in Yp c (to) has to be discarded or not. 
The situation for oxygen remains unclear since no oxygen 
abundances were measured at metallicities where the ef- 
fect is most pronounced ([Fe/H] < —3.0). However, in the 
range —3.0 < [Fe/H] < —1.5 there are many observations 
with [O/Fe] > 0.6 whereas the bulk of model stars in this 
range shows [O/Fe] ss 0.4. Furthermore, many observa- 
tions of Mg abundances in halo stars with [Fe/H] < —3.0 
exist, but only one shows a ratio of [Mg/Fe] > 1.0. The re- 
maining stars all have [Mg/Fe] < 0.8, which is in contrast 
to the predictions of the model. Contrary to O and Mg, 
there are indeed several observations of metal-poor halo 
stars with very high [Si/Fe] and [Ca/Fe] ratios at [Fe/H] 

< — 2.5 and the fit in [Si/Fe] and [Ca/Fe] is not too bad. 
However, there are some observed stars with [Si/Fe] < 0.0 
that are not reproduced by the inhomogeneous chemical 
evolution model. All told, model HI clearly does not fit 
the observations as well as model SI. 

H2: The iron yields shown in Fig. (model H2) are 
a result of the assumption that the 26 Mq model of SN 
1997D together with the 25 M Q model of SN 1998bw are 
correct. Coincidentally, Yp e (to) is also compatible with the 
10 M and 40 M Q models of SN 1997D and SN 1998bw 
due to the requirement to keep the average iron yield con- 
stant. Fig, [lq shows the resulting [el/Fe] distributions. Due 
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to the low amount of Ni ejected by the 26 Mq progenitor 
that is the same for models HI and H2 (c.f. Figs. O, H2 and 
Table |), the problems in [O/Fe] and [Mg/Fe] described 
in the discussion of model HI still persist. Compared to 
model HI, the fit in [Si/Fe] is now significantly improved, 
whereas model stars with [Ca/Fe] abundances that are 
too low are again generated by the inhomogeneous chem- 
ical evolution code (compare with model SI). However, 
although models HI and H2 predict metal-poor halo stars 
with [O/Fe] and [Mg/Fe] ratios as high as « 1.5, they 
can not be clearly discarded and the discovery of stars 
with metallicities [Fe/H] < -2.5 and 0.8 < [Mg/Fe] < 1.5 
would be a strong argument for the validity of the sudden 
decrease in the iron yields proposed by the models HI and 
H2. 

S2: Finally, in Fig. [l3| we show possible iron yields as- 
suming that the 10 M Q model of SN 1997D and 40 M 
model of SN 1998bw are correct (model S2). Here, Yp e (to) 
shows a plateau in the progenitor mass range from 15 Mq 
to 30 Mq , followed by an increasing yield with increas- 
ing progenitor mass. The resulting [el/Fe] distributions 
are shown in Fig. [17| A shallow slope in the [el/Fe] ratios 
is visible in every element, violating the constraint dis- 
cussed in Sect. (4.3.4). Furthermore, the scatter in [O/Fc] 
and [Si/Fe] clearly is not fitted by the model stars. This is 
especially conspicuous in the case of Si: According to the 
model we would expect the bulk of [Si/Fc] abundances 
to lie in the range 0.4 < [Si/Fe] < 0.8. To the contrary, 
most of the observed stars have [Si/Fe] < 0.4. Model S2 
therefore gives the worst fit to element abundance deter- 
minations in metal-poor halo stars. 

A physical explanation for a sudden drop of the 
iron yields of SNe with progenitor masse s around 25 Mq 
was suggested by Iwamoto et al. ( 2000 ): Observational 
and theoretical evidence indicate that stars with main- 
sequence masses M ms < 25 Mq form neutron stars with 
a typical iron yield of ~ 0.07 Mq, while progenitors more 
massive than this limit might form black holes and, due to 
the deep gravitational potential, have a very low (or no) 
iron yield. This might have been the case for SN 1997D. 
One of the two possible models reconstructing its light- 
curve assumes a 26 Mq progenitor and a very low kinetic 
energy of only a few times 10 50 erg and an equally low 
56 Ni yield of « 0.002 Mq. (But note that the lower- mass 
progenitor model seems more likely, c.f. Sect. 4.3.3). On 
the other hand, hypernovae such as 1997ef or 1998bw with 
progenitor masses around 30 Mq and 40 Mq and explo- 
sion energies as high as 10 — 100 x 10 51 erg might be 
energetic enough to allow for high iron yields even when a 
black hole forms during the SN event (see e.g. MacFadyen 
et al. 2001). However, one of the models of SN 1998bw 
proposes a 25 Mq progenitor with a kinetic energy typ- 
ical for hypernovae (i.e. much larger than the explosion 
energy of SN 1997D). This is in some sense a contradic- 
tion to the case of SN 1997D if we assume that a black 
hole formed in both cases: If the explosion mechanism is 
the same for SN 1998bw and SN 1997D it is natural to 
assume that the explosion energy scales with the mass of 



the progenitor and it is hard to imagine a mechanism that 
would account for a hypernova from a 25 Mq progenitor 
with an explosion energy that is w 100 times larger than 
the explosion energy from a 26 Mq progenitor. 

Therefore, model HI would fit nicely into the (qual- 
itative) hypernova scenario proposed by Iwamoto et al. 
( |2000| ), whereas models SI, S2 and H2 (yet) lack a physi- 
cal explanation. On the other hand, model SI gives a much 
better fit to the observations than HI, H2 and S2. Models 
HI and H2 could be tested, since they predict a number of 
stars with very high [O/Fe], [Mg/Fe], [Si/Fe] and [Ca/Fe] 
ratios (up to 1.5 dex) at metallicities [Fe/H] < —2.5. The 
discovery of such ultra a-element enhanced stars would 
be a strong argument in favour of the hypernova scenario 



proposed by Iwamoto et al. (2000) and the existence of a 
sudden drop in the iron yields of supernovae/hypernovae 
with progenitors around 25 Mq. 

Table || lists the numerical values of Yp e (to) as func- 
tion of progenitor mass to for the models discussed above. 
Model SI gives the best fit to the distribution of a-element 
abundances in metal-poor halo stars while S2 gives the 
worst. Although the models HI and H2 violate the con- 
straints discussed in Sect. ( 4.3.2; ) and thus do not give a 
fit as good as the one of SI, they cannot be ruled out on 
the basis of the observational data available to date. 



6. Conclusions 

Inhomogeneous chemical evolution models in conjunction 
with a current set of theoretical nucleosynthesis yields pre- 
dict the existence of very metal poor star s wit h subsolar 
[O/Fe] and [Mg/Fe] ratios (Argast et al. |2000| ). This re- 
sult is a direct consequence of the progenitor mass depen- 
dence of stellar yields, since core-collapse SNe of different 
masses imprint their unique clement abundance patterns 
on the surrounding ISM. No observational evidence of the 
existence of such stars is found, and recent investigations 
on the contrary indicate an increasing [O/Fe] ratio with 
decreasing metallicity [Fe/H] . This result of the inhomo- 
geneous chemical evolution calculations is primarily due 
to the input stellar yields and not due to the details of the 
model itself. This is a strong indication that the progen- 
itor mass dependence of existing nucleosynthesis models 
is not fully understood. This in itself is not surprising, 
since no self-consistent models of the core-collapse and 
the ensuing explosion exist to date (c.f. Liebendorfer et 



al. 2001; Mezzacappa et al. 2001; Rampp & Janka 2000). 
A crucial parameter of explosive nucleosynthesis models is 
the mass-cut, i.e. the dividing line between proto-neutron 
star and ejecta. This gives rise to a large uncertainty in 
the amount of iron that is expelled in the explosion of a 
massive star. On the other hand, oxygen and magnesium 
are mainly produced during hydrostatic burning and are 
therefore not strongly affected by the details of the ex- 
plosion mechanism. However, the distribution of [O/Mg] 
ratios of metal-poor halo stars suggests that either un- 
certainties exist even for O and Mg yields, or that obser- 
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Fig. 14. [el/Fe] distribution of model stars for O, Mg, Si and Ca resulting from the iron yields shown in Fig. [h] (model 
SI). 



vations overestimate oxygen or underestimate magnesium 
abundances in such metal-poor stars. 

The predictions of our inhomogeneous chemical evolu- 
tion model can be rectified under the assumption that the 
stellar yields of oxygen and magnesium reflect the true 
production in massive stars well enough, and by re plac- 
ing the stellar ir on yie lds of Thielemann et al. ([1996) and 
Nomoto et al. ( 1997 ) by ad hoc iron yields Yp e (m) as 
function of progenitor mass to. These are derived in this 
paper from observations of metal-poor halo stars and core- 
collapse SNe with known progenitor and ejected 56 Ni mass 



(the main source of 56 Fe by the decay 56 Ni — > 56 Co — » 
56 Fe). Such ad hoc iron yields have to satisfy the following 
constraints: First, the IMF averaged stellar yields should 
reproduce the mean [O/Fe] and [Mg/Fe] abundances of 
metal-poor halo stars. Second, the range and scatter ob- 
served in [O/Fe] and [Mg/Fe] ratios of metal-poor halo 
stars must be reproduced. This, in conjunction with stellar 
oxygen and magnesium yields, leads in turn to upper and 
lower boundaries for Yp c (m). Third, no slope should be 
introduced by Yp e (to) in the [el/Fe] distribution of model 
stars that is not compatible with observations. And finally, 
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Fig. 15. [el/Fe] distribution of model stars for O, Mg, Si and Ca resulting from the iron yields shown in Fig. [ll| (model 
HI). 



the progenitor mass dependence of the iron yields should 
be consistent with the ejected 56 Ni mass of observed core- 
collapse SNe with known main-sequence mass. Here, the 
situation is complicated by SN 1997D and SN 1998bw. 
The models recovering their light-curves give in each case 
two significantly different progenitor masses. These con- 
straints severely curtail the possible iron yield distribu- 
tions but are not stringent enough to determine Yp e (m) 
unambiguously. 

The main results of this paper are summarized in the 
following points: 



Observations of O and Mg abundances in metal-poor 
halo stars and of the ejected 56 Ni mass in core-collapse 
SNe, in conjunction with oxygen and magnesium yields 
from nucleosynthesis calculations and inhomogeneous 
chemical evolution models, provide a valuable tool to 
constrain the amount of iron ejected in a SN event as 
function of the main-sequence mass of its progenitor. 
The [el/Fe] distribution of model stars as function of 
metallicity [Fe/H] is sensitive to the iron yields of SNe 
with progenitors in the mass range 10—20 Mq. A steep 



increase of lp e (m) from w 1.5 • 10 Mq for a 10 M 
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Fig. 16. [el/Fe] distribution of model stars for O, Mg, Si and Ca resulting from the iron yields shown in Fig. [T^ (model 
H2). 



progenitor to w 5.5 • 10~ 2 Mq for a 15 Mq progenitor 
followed by a slow increase to « 7.0 • 10 _2 Mq for a 
20 Mq progenitor is required to give an acceptable fit 
to the observations of metal-poor halo stars. 
The further trend of Yp e (m) in the mass range 20 — 
50 Mq cannot be unambiguously determined by the 
available data. For this mass range we have deduced 
four possible iron yield distributions (models SI, S2, 
HI and H2) that explore the available freedom. These 
correspond to the four different combinations of prob- 
able progenitor masses of SN 1997D and SN 1998bw. 



Iron yield distributions that differ significantly form 
the presented models can be excluded. 
Model SI gives the best fit to observations while mod- 
els HI and H2 can not be ruled out. Model S2 gives 
the worst fit to observed [el/Fe] ratios in metal-poor 
halo stars. A change in the explosion mechanism of 
SNe II around 25 Mq is expected in the case of the 
"H" models. A test to distinguish between models SI 
and H1/H2 would be the discovery of very metal-poor 
stars ([Fe/H] < —2.5) that are highly enriched in a- 
elements. 
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Fig. 17. [el/Fe] distribution of model stars for O, Mg, Si and Ca resulting from the iron yields shown in Fig. \YQ (model 
S2). 



5. Iron yield distributions derived from observations 
through inhomogcncous chemical evolution models 
yield constraints on the mass-cut in a SN II event 
if the detailed structure of the progenitor model is 
known (i.e. the size of the iron core and the zone 
that undergoes explosive Si burning). Thus, they can 
be used as benchmarks for future core-collapse super- 
nova/hypernova models. 

In future, a large and above all homogeneously ana- 
lyzed sample of O, Mg and Fe abundances in very metal- 



poor stars is needed to derive more stringent constraints 
on Ype (tti) . Not only would this allow us to determine the 
exact extent of the scatter in [O/Fe] and [Mg/Fe] ratios, 
but would also answer the important question whether the 
scatter in [O/Mg] is real or due to some (yet) unknown 
systematic errors in O and Mg abundance determinations. 
If the scatter in [O/Mg] turns out to be real, updated 
nucleosynthesis calculations including rotation and mass 
loss due to stellar winds are needed to understand O and 
Mg abundances in metal-poor halo stars. 
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Also very valuable would be the observation and anal- 
ysis of further core-collapse supernovae/hypernovae. Only 
six core-collapse SNe with known progenitor and ejected 
56 Ni mass are known to date, and for two of them their 
progenitor masses are not clearly determined. Especially 
the discovery of a SN II with a progenitor in the critical 
mass range from 20 — 30 Mq could provide us with the 
information needed to discern between the four models 
presented above, or at least whether the "S" or "H" mod- 
els have to be preferred. This would also be a step towards 
answering the question whether a change in the explosion 
mechanism of core-collapse SNe occurs, i.e. the formation 
of a black hole and significant increase of the explosion 
energy for m > 25 Mq . 
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